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ABSTRACT 

We study the embedding of D(8-p)-branes in the background geometry of parallel 
Dp-branes for p < 6. The D(8-p)-brane is extended along the directions orthogonal to the 
Dp-branes of the background. The D(8-p)-brane configuration is determined by its Dirac- 
Born-Infeld plus Wess-Zumino-Witten action. We find a BPS condition which solves the 
equation of motion. The analytical solution of the BPS equation for the near-horizon and 
asymptotically flat geometries is given. The embeddings we obtain represent branes joined 
by tubes. By analyzing the energy of these tubes we conclude that they can be regarded as 
bundles of fundamental strings. Our solution provides an explicit realization of the Hanany- 
Witten effect. When p = 6 the solution of the BPS equation must be considered separately 
and, in general, the embeddings of the D2-branes do not admit the same interpretation as 
in the p < 6 case. 
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1 Introduction 



D-branes have been shown to play a fundamental role in the non-perturbative structure of 
string theory. Originally [|I|, they were introduced as hypersurfaces in space-time on which 
open strings with Dirichlet boundary conditions are allowed to end. The extended objects so 
defined are charged under the Ramond-Ramond (RR) sector of the theory and are essential 
to understand the duality structure of string theories. Moreover, the low energy physics of 
D-branes can be described by supersymmetric gauge theory and, actually, super Yang-Mills 
theory arises as a low energy description of parallel D-branes. 

D-branes also appear as solutions of the supergravity effective action of string theory [0. 
These solutions are extended generalizations of the Schwarzschild geometry which describe 
the gravitational field created by an object carrying Ramond-Ramond charge. This gravita- 
tional field acts on the dimensions orthogonal to the D-brane worldvolume which, from the 
point of view of the gauge theory, can be regarded as "extra dimensions" . 

The interplay between gauge theory and gravity is on the basis of the Maldacena conjec- 
ture 0. Indeed, Maldacena has argued that there is a remarkable duality between classical 
supergravity in the near-horizon region and the large N 't Hooft limit of the SU(N) Yang- 
Mills theory. In particular, the near-horizon geometry of parallel D3-branes is equivalent to 
the space AdS§ x S 5 , whose boundary can be identified with four-dimensional Minkowsky 
space-time. In this case one gets the so-called AdS/CFT correspondence between type IIB 
superstrings on AdS§ x S 5 and M = 4 non-abelian super Yang-Mills theory in four dimensions 

[5], [6|]. Within this context, one can compute Wilson loop expectation values by consider- 
ing a fundamental open string placed on the interior of the AdS§ space and having its ends 
on the boundary J7|. This formalism has been used to extract quark potentials both in the 
supersymmetric and non-supersymmetric theory. Moreover, Witten || has proposed a way 
to incorporate baryons by means of a D5-brane wrapped on the S 5 , from which fundamental 
strings that end on the D3-branes emanate || [UJ. 



In ref. ||1 1|| , the problem of a D5-brane moving under the influence of a D3-brane back- 
ground was considered. The D5-brane dynamics is governed by an action which is the sum 
of a Dirac-Born-Infeld and a Wess-Zumino term. The equations of motion of the static 
D5-brane are solved if one assumes that the fields satisfy a certain first-order BPS differ- 
ential equation that was found previously by Imamura Hl2| . This BPS condition can be 



integrated exactly in the near-horizon region of the background D3-brane geometry and the 
solutions are spikes of the D5-brane worldvolume which can be interpreted as bundles of 
fundamental strings ending on the D3-branes. Actually, these kinds of spikes are a general 
characteristic of the Dirac-Born-Infeld non-linear gauge theory and can be used to describe 



strings attached to branes [13, 14 



The BPS condition for the full asymptotically flat D3-brane metric was analyzed nu- 



merically in ref. [11]. As a result of this study one gets a precise description of the string 



creation process which takes place when two D-branes cross each other, i.e. of the so-called 



Hanany- Witten effect [15, [lq| . These results were generalized in ref. [17| to the case of 



baryons in confining gauge theories. Moreover, in ref. |Tj| the BPS differential equation was 



obtained as a condition which must be fulfilled in order to saturate a lower bound on the 
energy. 
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In this paper we shall study the worldvolume dynamics of a D(8-p)-brane moving un- 
der the action of the gravitational and Ramond-Ramond fields of a stack of background 
Dp-branes for p < 6. The D(8-p)-brane is extended along the directions orthogonal to the 
worldvolume of the background Dp-branes. For p = 3 this system is the one analyzed in ref. 



IIJ . By using the method of ref. ||18|| , we will find a BPS first-order differential equation 
whose solutions also verify the equation of motion. Remarkably, we will be able to integrate 
analytically this BPS equation both in the near-horizon and asymptotically flat metrics. 



For p < 5 the solutions we will find are similar to the ones described in ref. JTTJ, i.e. the 



D(8-p)-brane worldvolume has spikes which can be interpreted as flux tubes made up of 
fundamental strings. Our analytical solution will allow us to give a detailed description of 
the shape of these tubes, their energy and of the process of string creation from D-brane 
worldvolumes. 

The organization of this paper is the following. In section 2 we formulate our problem. 
The BPS differential equation is obtained in section 3. In section 4 this BPS condition is 
integrated for p < 5 in the near-horizon region. The p = 6 system requires a special treat- 
ment, which is discussed in section 5. The solution of the BPS equation for the full Dp-brane 
geometry is obtained in section 6. The detailed analysis of this solution is performed in ap- 
pendix A. Finally, in section 7, we summarize our results and point out some directions for 
future work. 



2 Worldvolume Dynamics 



The ten-dimensional metric (in the string frame) corresponding to a stack of N coincident 
extremal Dp-branes in type II superstring theory is H: 



ds 2 



H P (r) 



-dt 2 + dxl) + 



H p (r) 



( dr 2 + r 2 dSll_ p ) . 



(2.1) 



In eq. ( |2.1| ) x \\ are p-dimensional cartesian coordinates along the branes, r is a radial coordi 

I 2 

L 8-p 



nate parametrizing the distance to the branes and dQl is the line element of an unit 8 — p 



sphere. The harmonic function H p (r) appearing in eq. ( [2.1|) is given by: 

R\7- P 



H P {r) 



a + 



(2.2) 



where a = 0, 1 and the 'radius' R can be written in terms of the number N of branes, the 
Regge slope a ' and the string coupling constant g s as follows: 



R 7 ~ p 



(2.3) 



We have chosen our conventions in such a way that the tension of the fundamental string Tf 
is: 



T 



1 



2-na' 



(2.4) 



The value a = 1 in eq. (|2.2|) corresponds to the full geometry of the stack of Dp-branes. 
Actually, in order to have an asymptotically flat space-time, we shall restrict ourselves to 
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the case p < 7. Taking a = in the harmonic function H p (r) is equivalent to approximate 



the Dp-brane geometry by the near- horizon metric of the 'throat' region of (2.1). 

The Dp-brane solution we are considering is also characterized by a dilaton field <j>(r) and 
a Ramond-Ramond (RR) p + 2-form field G. The corresponding values are: 



(r) 



H P (r) 



P-3 
4 



G = *- 

'"' ,x \\ ,r dr 



H P (r) 



(2.5) 



where 0(r) = 0(r) — 0^ and 0oo is the value of the dilaton field at infinite distance of 
the Dp-branes (i.e. 0oo = lim r ^oo <f>(r)). As was pointed out in ref. ||, in order to trust 
this type II supergravity solution, both the curvature in string units and the dilaton must 
be small. This fact introduces restrictions in the values of the radial coordinates for which 
the correspondence between the supergravity and gauge theory descriptions is valid (see also 
ref. [0). 

Let us now consider a D(8-p)-brane embedded in the transverse directions of the stack of 
the background Dp-branes. The dynamics of this D(8-p)-brane is determined by its action, 
which is a sum of a Dirac-Born-Infeld and a Wess-Zumino term: 

S = -T 8 _ p J d 9 ~ p i e~$ y/-det(g + F) + T 8 _ p / d 9 ~ p £ AA*G, (2.6) 

where g is the induced metric on the worldvolume of the D(8-p)-brane, A is a worldvolume 
abelian gauge field and F = dA its field strength. In eq. (|2.6|) *G is a (8-p)-form, which is 
the pullback of the Hodge dual of the background RR (p+2)-form G. Notice that this RR 
field acts as a source for the worldvolume gauge field A. The coefficient Tg_ p appearing in 
the action fl2.6|) is the tension of a D(8-p)-brane, which is given by 0: 

T 8 _ p = (2ny~ 8 (a') E ^(g s )- 1 . (2.7) 

Let 9 1 , 8 2 , ■ ■ 8 8 ~ p be coordinates which parametrize the S s ~ p transverse sphere. The 
worldvolume coordinates £ a (a = , • • • , 8 — p ) will be taken as: 

r = (t,6\--- ,9*-P) ■ (2-8) 

We shall assume that the #'s are spherical angles on S s ~ p and that 9 = 9 S ~ P is the polar 
angle (0 < 9 < n). Therefore, the S 8 ~ p line element dVl\_ p can be decomposed as: 

dn 2 8 _ p = d9 2 + (sinfl) 2 dQ 2 7 _ p , (2.9) 

where dfl 2 _ p , which only depends on 9 1 , ■ ■■ ,9 7 ~ p , is the line element of a S 7 ~ p sphere. 
We are going to consider static configurations of the D(8-p)-brane which only depend on the 
polar angle 9. Actually, we shall restrict ourselves to the case in which the only non-vanishing 
fields are r(9), x(9) and Aq(9), x being a direction parallel to the background Dp-branes. It 
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is straightforward to verify that the action for such configurations can be written as: 
S = TT 8 _ p fi 7 _ p / dS (sin0) 7 ~ p x 



x 



-r 7 - p H p (r) 



\ 



■ 12 



ry 2 j ry f 2 j 



H P (r) 



- Flo + ("1) P+1 (7 - p) A R 7 ~ p 



(2.10) 



where the prime denotes derivative with respect to 9 and T = J dt. In eq. (|2.10|) fi 7 -p is 
the volume of the unit (7-p)-sphere, given by: 



n 



7-P 



2vr 2 



(2-11) 



It is interesting at this point to remember that we are considering static solutions of the 
gauge field that only depend on 9. The electric field for these configurations is: 



E = F 0t8 = -de A . 
Let us now define the displacement field D p as: 



D p {9\ 



-l) p+1 OS _ ^ +1 (sm9) 7 ~ p r 7 - p H p (r)E 



T T 8 _p f2 7 _ p dE 



ry'2i j ^*^2 j 



H P {r) 



(2.12) 



(2.13) 



The extra factors appearing in the definition ( 2.13Q have been included for convenience. The 
Euler-Lagrange equation for Aq implies an equation for D p which is easy to find, namely: 



d 
d9 



D p (9) = -(7-p)R 



sin 



0)7 -p 



(2.14) 



Notice that the right-hand side of eq. ( |2.14j ) only depends on 9. Therefore ( |2.14j ) can be 
integrated and, as a result, one can obtain D p as a function of 9. We shall work out the 
explicit form of D p {9) at the end of this section. At present we only need to use the fact 
that D p ( 9 ) satisfies eq. (gig) . Actually, substituting ( 7 - p ) R 7 ~ p ( sin 9 ) 7 ~ p by -dg D p in 
the Wess-Zumino term of the action and integrating by parts, one can recast S as: 



S 



-TU 



(2.15) 



with U given by: 



U 



d9 



-1) P+1 ED P (9) + 



+ (sin#) 7 - p r 7 - p 



H P (r) 



\ 



■ 12 



/y 2 j ry I '2 j 



H p {r) 



E 2 



. (2.16) 
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Notice that, as is evident from their relation ( 2.15 ), S and U give rise to the same Euler- 
Lagrange equation. Actually, since we have eliminated A in favor of D p , we can regard U 
as the Legendre transform of S and, therefore, U can be considered as an energy functional 
for the embedding of the D(8-p)-brane in the Dp-brane background. The fields E and D p 
are related, as is obvious from eq. ( |2.13| ). It is not difficult to invert eq. (|2.13| ) and get E in 



terms of D p . The result is: 



E 



-iy +l 



iy 2 j iy f 2 j 



H P {r) 

\ (D p (9)) 2 + [(sm6y-Pr 7 ~PH p 



r 



D p {9) . 



(2.17) 



Using this relation we can eliminate E from the expression of U: 



U 



p ^ L 7—p 



■ 12 



+ 



X 



12 



H P {r) 



(D p (8)) 2 + [(sm 9) 7 -Pr 7 'PH p (r) 



(2.18) 



Recall that D p (8) is a known function of 9 (see below). Therefore eq. fl2.18|) gives the 
energy functional U in terms of x(8) and r(6). These functions must be solutions of the 
Euler-Lagrange equations obtained from U. From the study of the functional U in several 
situations we will determine the shape of the D(8-p)-brane embedding in the background 
geometry. Let us consider, first of all, the near-horizon approximation, which is equivalent 



to taking a 
given by: 



in the harmonic function ( |2.2| ). In this case r 7 p H p (r) = R 7 p and U is 



U 



T 8 „ P 1V P fdeJr 2 + r' 2 + 



r 7-p 

R 7 -p 



x 



12 



(D p (6)) 2 + (smOfV-P) RW-p) . (2.19) 



By inspecting eq. Q2.19] ) one easily concludes that U transforms homogeneously under a 
simultaneous rescaling of the radial and parallel coordinates of the form: 



[r,x] 



ar, a 2 x 



(2.20) 



where a is a constant. It follows that the equations of motion derived from U are invariant 
under the transformation ( |2.20| ). Moreover, the homogeneous character of U under the 
transformation (|2.20|) implies the following scaling law: 



x 



R 2 

5— p 5 
f 2 



(2.21) 



where the power of R has been determined by imposing dimensional homogeneity of both 
sides of fl2.21|) . Eq. (|2.21| ) is precisely the holographic UV/IR relation found in ref. [20|. 
According to eq. Q2.21| ), for p < 5, large radial distances (r — > oo) correspond to small 
values of the parallel coordinate x. For p — 5, the distance x is insensitive to changes in 
r, whereas, for p = 6, x increases when r grows. The consequences of this relation for the 
correspondence between field theories and near-horizon supergravities have been discussed 
in ref. [TJUJ (see also ref. PIJ). In our approach we shall see that, indeed, the p = 6 case is 
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special and the embedding of the D2-brane in the D6-brane background geometry has new 
characteristics, which must be studied separately. 

Let us finish this section by giving the expressions of the displacement fields D p . The 
^-dependence of D p can be obtained by integrating the right-hand side of eq. (|2.14j ). The 
results one gets for the different cases are: 



D {9 

D 2 (6 

D 3 {6 
D 4 (6 
D 5 (9 
D 6 (9 



R 7 



R b 



R 5 



R 4 



cosfl(sin 6 # + -sin 4 # + -sin 2 # + — ) + — (2u - 1 

5 5 5 5 



5 15 15 

cos 9 ( sin 5 6> + - sin 3 ^ H sin 9 ) H (ttv - 9 

v 4 8 ; 8 v 



cos 9 ( sin 4 fl + - sm 2 9 + -) + -( 2v - 1 
3 3 3 



cos 9 ( sin 3 9 + - sin 9 ) + - ( nu - 9 



(2.22) 



iT cos 9 ( sin 2 fl + 2 ) + 2 ( 2u - 1 



R 2 cos 9 sin 9 + -nu — 9 



R 



cos 9 



2u - 1 



where we have parametrized the additive constant of integration by means of a parameter 
v. We have chosen the ^-independent term in D p {9) in such a way that the value of the 
displacement field at 9 = n is given by: 



D p (n) = -2v^F )J { R 7 ~ p (l 



r(V) 



(2.23) 



Eq. ( gjg) allows to give a precise meaning to the parameter v |[L1|| . In fact, as we shall 
verify below, there exist solutions for which r — ► oo and x' — > when 9 — > 7r in a way that 
simulates a 'flux tube' attached to the D(8-p)-brane. The 'tension' (i.e. the energy per unit 
radial distance) for one of these spikes can be obtained from the expression of U in eq. ( |2.18| ). 
Actually, as 9 — > tt one can check that the terms with r ' dominate over the other terms in 
the first square root in eq. ( |2.18[ ), whereas the second square root can be approximated by 
| D p (tt) |. As a result ||11|| , the tension of the spike is Ts- P fl7- P \ D p (tt) |. Using the value of 
D p (n) given in eq. (|2~23|) and the fact that (see eqs. Q, Q and 



NT f r( 7 -^ 



(2.24) 



we conclude that the tension of the 9 = tt spike is (1 — v) NTf, where 7/ is the tension of 
the fundamental string written in eq. (|2.4|). This result implies that we can interpret the 
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9 = 7r tube as a bundle of (l — u)N fundamental strings which, from the gauge theory point 
of view, corresponds to a baryon formed by (1 — v) N quarks. It is clear that, although in 
the classical theory v is a continuous parameter, upon quantization v should be a multiple 
of 1/N taking values in the range < v < 1. 

Let us finally point out that there exist other solutions for which r — > oo and x' — > 
when 9 —>■ 0. The asymptotic tension for these solutions is T 8 _ P Q 7 _ P | D p (0) |. The value of 
D p (9) at = can be obtained from the values given in eq. ( |2.22j ). The result is: 

DM = 2V^^-^R 7 - p v. (2.25) 



Using again eq. (|2.24|) , one can verify that T S _ P Q 7 _ P \ D p (0) \ = v NTf, which corresponds 



to a bundle of v N fundamental strings. This fact provides an interpretation of v for this 
second class of solutions. 



3 BPS conditions 

In the remaining of this paper we are going to study solutions of the equations of motion of 
the brane probes which are not extended in the directions parallel to the background branes. 
This amounts to take x' — (i.e. x(9) = constant) in our previous equations. From the 
expression ( |2.18| ) of U we can get the differential equation which determines r as a function 



of 9. Indeed, the Euler-Lagrange equation derived from U is: 



d 
d9 



.^J ty 2 j ty ~t 2 



(D p (9)) 2 + [(sm9) 7 -Pr 7 -PH p { 



^ [D p {9)f + [(sm9y-Pr^H p (r)f + (3.1) 

a / ty * _L ty I & * 

Vr 2 + r' 2 (r sin 9 H p (r) 



+ (7-p)a 



{D p (9)) 2 + [{sin 9) 7 -Pr 7 -PH p {r) 



Trying to obtain a solution for such a complicate second-order differential equation seems, 
a priori, hopeless. A possible strategy to solve eq. fl3.1| ) consists of finding a first integral 
for this system. Notice that the integrand of U depends explicitly on 9 (see eq. ( |2.18| )) 
and, therefore, there is no first integral associated to the invariance under shifts of 9 by a 
constant. Nevertheless, we will be able to find a first-order equation such that any function 
r(9) satisfying it is a solution of the equation ( |3.1|) . This first-order condition is much simpler 
than the Euler-Lagrange equation (|3.1|) and, indeed, we will be able to solve it analytically, 
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both in the near-horizon (a = 0) and asymptotically flat (a = 1) geometries. For p = 3 and 



a = 0, the first order equation was found by Imamura [|12j as a BPS |22fl condition for the 
implementation of supersymmetry in the worldvolume theory of a D5-brane propagating in 
a D3 background. This condition was subsequently extended to (p = 3, a — 1) and (p = 4, 
a = 0) in refs. JTT[] and []T7| respectively. In order to generalize these results, we shall follow 
here the approach of ref. 



23J), where the BPS equation was found by 
tb 

let us define the quantity: 



|T8| (see also ref 

requiring the saturation of a certain bound on the energy functional 
Following ref. 



A p (r. 



9 



r 7 - p H p (r) 



sin 9 



In terms of A p ( r , 9), we define the function f p (r, 9) as follows: 

- D p (6)cos6 



fp(r,e) 



A p (r, 6) sm6 
A p (r , 9) cos 9 



D p (9)sin9 ' 



(3.2) 



(3.3) 



Making use of the function f p it is easy to verify that the energy can be put in terms of a 
square root of a sum of squares. Indeed, it is a simple exercise to demonstrate that U can 
be written as: 



U 



p ° L 7—p 



d9 x 



x yz 2 + r 2 

where Z is given by: 



A p (r, 9)cos9- D p (9)s\n9\ [--f p {r,\ 



(3.4) 



Z = r 



A p (r, 6) cos 9- D p (6) sin0 l + -/ p (r,^) 



In view of eq. 



it is clear that the energy of the D(8-p)-brane is bounded as: 



U > T» r ,n 



8-p "7-p 



d9 



Z 



(3.5) 



(3.6) 



This bound is saturated when r(9) satisfies the following first-order differential equation: 



f P (r, 9) . 



(3.7) 



It is straightforward to prove that any function r(9) that satisfies eq. (|3.7| ) is also a solution of 
the Euler-Lagrange equation (|3.1|) . Notice that the condition ( |3.7p involves the displacement 
field D p {9) (see eq. (|3.3| )). The expressions of the D p s have been given at the end of section 
2 (eq. ( |2.22| )). However, in order to demonstrate that the solutions of eq. ( |3.7| ) verify the 
equation of motion ( |3.1|) , we do not need to use the explicit form of D p {9). The only property 
required for this proof is the value of the derivative of D p (9), given in eq. ( [2.14| ). Moreover, 
as pointed out in ref. |18[ for the p = 3 case, Z can be written as a total derivative: 



8 



Z = Te l D v^ rcosd + (g^ + ^T) ("i^rl • (3.8) 

It is important to stress the fact that eq. (|3.8|) can be proved without using the condition 
( j3.7| ) or the Euler-Lagrange equation ( |3.1| ) (only eq. ( j2.14|) has to be used). Eq. ( |3.8|) implies 
that the bound (|3.6|) does not depend on the detailed form of the function r(9). Actually, as 
a consequence of eq. ( |3.8| ), only the boundary values of r(9) matter when one evaluates the 
right-hand side of eq. (|3.6| ). The functions r(9) which solve eq. (|3.7|) correspond to those 
D(8-p)-brane embeddings which, for given boundary conditions, have minimal energy. Due 
to this fact we shall refer to eq. (|3.7|) as the BPS condition. The search for the solutions of 
the BPS equation and their interpretation will be the subject of the next three sections. 

4 The near-horizon solution for p < 5 

In this section we are going to solve the BPS first-order differential equation in the near- 
horizon geometry for p < 5. It will become clear in the process of finding the solution that 
the p = 6 case is singled out (as expected from the holographic relation fl2.21| )). This p = 6 
case will be discussed separately in section 5. 

Let us start by defining the function A p ( 9 ) by means of the equation: 



AJ9) = — — 
pV J RT - P 



R 7 - p (sin9) 6 ~ p cos9 - D p (9 



(4.1) 



The right-hand side of eq. (|4.1| ) is known and, thus, A p ( 9 ) is a known function of 9 whose 
explicit expression can be obtained by substituting the values of D p { 9 ), given in eq. ( |2.22| ), 
in (|4.1|) . Moreover, the BPS condition ( p.7|) for a = can be written in terms of A p ( 9 ) as: 

r' (sin#) 6 -P - A p (#)cos# 

7 ~ A p (#)sin# ' ^ ' ' 

A key point in what follows is that A p (9) has a simple derivative, which can be easily 
obtained from eq. (|2.14j) , namely: 

±A P (9) = (6-p)(sm9) 5 ^. (4.3) 

By inspecting the right-hand side of this equation we observe that p = 6 is a special case. 
Indeed, when p ^ 6, we can represent the (sin#) 6-p term appearing in the BPS condition 
as: 

sin 9 n 

(sin#) 6 - p = — W K(9), (p^6). (4.4) 

After doing this, the right-hand side of eq. (|4.2|) is immediately recognized as a logarithmic 
derivative and, as a consequence, the BPS condition can be readily integrated. The result 
one arrives at is: 
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Figure 1: The functions A p (9) and D p (6) for < 9 < n. For illustrative purposes we have 
plotted these functions for p = 4 and v = 1/4. 



C 



sin 9 



i 

6— p 



(4.5) 



where C is a positive constant. For p < 5 there is a fractional power of A p ( 9 ) in the right- 



hand side of eq. ( flop and, thus, the solution we have found only makes sense for those values 
of 9 such that A p ( 9 ) > 0. Moreover, it is clear from eq. ( [4.3| ) that: 



<iA p 



> 



> < 6, < # < 7T 



(4.6) 



and, therefore, A p (9) is a monotonically increasing function in the interval < 9 < n. The 



values of A p ( 9 ) at 9 



0, 7r are: 



A p (0) 



A p (tt) 



A»(Q) 

R 7 -p 



DM 

R7-p 



2v^ ) 7 2 ( ^ < , 



1 - u) > 



(4.7) 



where, in order to establish the last inequalities, we have used the fact that < v < 1. From 
this discussion it follows that there exists a unique value 9q of the polar angle such that: 



A p (# ) = 0. (4.8) 

The form of the function A p (9) has been displayed in figure 1. It is clear that the solution 
( [1.5D is valid for 9 < 9 < n. Notice that 9 depends on v and, actually, is a monotonically 
increasing function of v. In fact, from eq. ( |4.7|) one concludes that 9 = for v = 0, 
whereas 9q = tt for u — 1. 



For p = 5 the function appearing in the right-hand side of eq. ( f4.5|) is: 

A 5 (9) = 9 - . 



(4.9) 
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Thus, our solution in this case is: 



9 



C 



IT V 



sin 9 



(4.10) 



As the radial coordinate r must be non-negative, from eq. ([4.1 0|) one immediately concludes 
that, also in this p = 5 case, 9$ < 9 < ir, where now 9q depends linearly on z/, namely: 



9n 



7T V 



V 



(4.11) 



The solution we have found coincides with the one obtained in refs. [IT, |T7]] for p — 3, 4. Our 
result (|4.5|) generalizes these solutions for any p < 6 (the solution for p = 6 will be given in 
the next section). By inspecting eq. ([15]) it is easy to conclude that, for v ^ 1, riff) diverges 
when 9 — > n. Actually, when 9 ?a n and v ^ 1, r(#) behaves as: 



r{9) 



C 



A p (vr) 



1 

6—p 



7T 







(1/ ^ 1 



(4.12) 



which corresponds to a "tube" of radius C A p ( tt ) 6 p . We will check below that the energy 
of these tubes corresponds to a bundle of (1 — v)N fundamental strings emanating from the 
D(8-p)-brane. When v ^ the critical angle 9q is greater than zero and, as a consequence, 
eq. (|4.5|) gives r(6*o) = 0. Actually, when v ^ 0, the solution has the shape represented in 
figure 2a. On the contrary, when v = 0, the function r(9) takes a non- vanishing value at the 
critical value #o = 0. Notice that, in this last case, both numerator and denominator on the 
right-hand side of eq. ( |4.5| ) vanish. The behaviour of A p (#) near 9 = can be obtained by 
a Taylor expansion. The successive derivatives of A p (9) can be easily computed from the 
value of its first derivative, given in eq. ( f4.3| ). It is easy to prove that the first non- vanishing 
derivative of A p ( 9 ) at 9 = is: 



d9 e ~P 



AJ9 



6 - p] 



8=0 



Moreover, the (7 — p) derivative of A p ( 9 ) at 9 = vanishes, i.e. 



. 



(4.13) 



(4.14) 



6»=0 



Therefore, the expression of A p ( 9 ) for values of 9 close to zero takes the form: 



A p (9) w A p (0) + 9 b ~ p + c p 9*~ p + ■•■ , (4.15) 

where c p is a non- vanishing coefficient. Substituting this expansion in the right-hand side of 
eq. ( |4.5| ) for v = leads to the conclusion that: 

r(0) = C, (* = 0), (4.16) 

ie. r(0) 7^ for z/ = 0, as claimed above. Moreover, r'( ) vanishes in this case: 

r'(0) = 0, (i/ = 0) . (4.17) 
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I 
I 



I 
] 



(a) 



(b) 



Figure 2: Plot of the near-horizon solution Q4.5|) for (p = 3, v = 1/4) (a) and (p = 3, v = 0) 
(b). The discontinuous line represents the polar axis {9 = n at the top). The dot corresponds 
to the origin r = 0. The background branes are located at the origin and extend in the 
directions orthogonal to the plane of the figure. 



The shape of the solution for v = has been plotted in figure 2b. Notice that the v = 
tube corresponds to the ordinary baryon with iV quarks. 

The solution of the near-horizon BPS condition written in eq. ( |4.5|) is not the only 
one. Following ref. |I1J, we can construct a new solution, valid for values of 9 in the range 
< 9 < 9q, in which AJ 9 ) < 0, as follows: 



9 



C 



-A p (6) 



6—p 



sin 9 



(4.18) 



Eq. ([4.1 8|) describes a spike at 9 — 0. We shall refer to the solution ( |4.18| ) as a "lower tube" 
solution, in contrast to the "upper tube" solution of eq. ( |4.5|) . Actually, these two types of 
solutions are related. In order to verify this fact, let us point out that the functions D p (9) 
and A p ( 9 ) change their sign under the transformation 6 1 — > 7r — 9, v — > 1 — v: 



D p (9; v) 



-D p (n -9; 1 - v) , 



A p (9; v\ 



-A p (tt -9; 1 -v 



and, therefore, the solutions fl4.5|) and ( |4.18|) are simply related, namely: 



(4.19) 



f (9; v) = r(7r - 9; 1 - v) . 
It follows from this relation that the lower flux tubes correspond to vN quarks. 



(4.20) 
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Let us now evaluate, following ref. [18|, the energy for the two types of near- horizon 
solutions we have found. These solutions saturate the bound ( |3.6|) and, thus, the energy is 
precisely the right-hand side of this equation. Moreover, when r is a solution of eq. ( |3.7| ), Z 
can be written as: 



Z = r[A p (r, 0) cos9 - D p (6)sm6][l + (J—) 2 



(4.21) 



It is obvious from this equation that, for our solutions, the sign of Z is just the sign of 
A p (r, 9) cos 9 — D p (9) sin 6*. Moreover, using the definitions of A p (eq. ( |3~2]) ) and A p 
(eq. (|4.1|)), one can prove that: 



A p ( r , 9) cos ( 



£> p (0)sin0 



ar 7 p (sin6>) 7 p cos( 



R J - p A p (9)sm9 



(4.22) 



where we have used the general harmonic function fl2.2p . In the near-horizon case (a = 0) 
the first term on the right-hand side of eq. ( |4.22|) is absent and, therefore, the sign of 
A p (r , 9 ) cos6> — D p ( 9)sin9 is just the sign of A p (9). In the solutions ( pT5| ) and fl4.1cf ), the 
angle 9 can take values in a range such that A p (9) has a well defined sign. Therefore, we 
can write: 



d9 



Z 



d9Z 



If we now define Z as: 



Z = T 8 _ p r27_ p J d9 Z , 
it is immediate that the energy of our near-horizon BPS solutions is: 



(4.23) 
(4.24) 



Ubps = \Z\. (4.25) 

In order to evaluate Z for the near-horizon solutions, we make use of the representation 
( P-8p of Z as a total derivative. Recall that 9 varies in the range 9i < 9 < 9 f where 9i = 9 
{9i = 0) and 9f = tt (9f = 9o) for the upper (lower) tube solution. It follows that Z can be 
written as the sum []: 

Z = Z s + Z gs , (4.26) 

where Z s and Z gs are given by: 



T 8 - p fi 7 - P D p {9)r{9) cos 9 



T s . p n 7 . p R 7 ~ p r{9) 



sin^ 



(4.27) 



It is clear that Z only depends on the values of r(9) at the boundaries 9 = 9{ and 9 = 9f. In 
this sense we can regard Z as a topological quantity which, for fixed boundary conditions, 
is invariant under local variations of the fields 



1 The labels s and gs refer to "spike" and "ground state", according to the interpretation given to these 
two contributions in ref. |18l. 
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Let us now compute Z ga . After substituting the solutions 
equations, we get: 



and (|OD of the BPS 



Ts-p R 7 p C 



sm i 



\7-p 



±A„(# 



1 

6— p 



(4.28) 



where the + (— ) sign corresponds to the upper (lower) tube solution. The angles 9{ and Of 
can take the values 0, 7r and 6*0. For = 0, n the right-hand side of eq. (|4.28|) vanishes due to 
the ( sin i 
(see eq. 



) 7 p factor, whereas = gives also a vanishing contribution because A p ( ) =0 
8h). In conclusion, we can write: 



'gs 







(4.29) 



In the same way one can obtain the values of Z s . One can check in this case that the 
contribution of = 0q to the right-hand side of the first equation in (|4.27|) vanishes and, as 
a consequence, only 6 = 0, 7r contribute. After using eqs. (|2.23| )- (|2.25| ), one gets the result: 



T 



NTf L , (upper tube) 



vNTf L 



(lower tube) 



(4.30) 



where L = r(n) (L = r(0)) for the upper (lower) tube solution. It is now evident that 
Ubps — \ Z s \ is equal to the energy of a bundle of (1 — u)N or uN fundamental strings. 
Notice that this result is the same as the one found at the end of section 2 for the energy of 
the 6 = and 9 = n spikes. Let us finally mention that in ref. JTBJ it has been argued in 
favor of interpreting Z s as a central charge in the worldvolume supersymmetry algebra. 



5 The near-horizon D6-D2 system 

In this section we are going to integrate the BPS condition for p = 6. Notice that, according 
to the definition (|4.1|) and the expression of Dq{6) given in eq. Q2.22j ), A 6 (9) is constant, 



i.e.: 



Ah(9' 



1 - 2u . 



(5.1) 



For p = 6 we cannot use eq. 



and, hence, we have to deal directly with eq. 



Actually, eq. (|4.2f) makes sense only for v ^ 1/2 since Aq ( 9 ), and therefore the denominator 



of eq. 
is: 



2D vanishes identically for v = |. With this restrictionF], the p = 6 BPS condition 



r 
r 



1 



1 



cos 9 



1 — 2v sin# sin# 

It is not difficult to realize that the right-hand side of eq. 
derivative, namely: 



(5.2) 

721) can be written total 



2 The fact that the p 
approximation. 



6, v = 1/2 BPS condition is ill-defined is an artifact of the near-horizon 
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r 

T 



d 



log 



(tan § )' 



sin 9 



(5.3) 



Therefore, the BPS equation can be immediately integrated. It is convenient to put the 
result in the form: 



r(9) = A 





sin f 


1 - 2u 







cos • 



(5.4) 



where A is a positive constant. The nature of the solution Q5.4D depends critically on the 
value of v. First of all, the range of values of 9 is not restricted, i.e. < 9 < n. If < v < 1/2, 
the cos (9/2) term present in the denominator of (|5.4| ) makes r — > oo at 9 w 7r, whereas for 
1/2 < i/ < 1 the solution diverges at # ~ 0. Actually, the solution for < v < 1/2 is related 
to the one for l/2<z/<lby means of the equation: 



r(9; u) = r(n - 9; 1 - u) 



(5.5) 

Eq. ( |5.5| ) can be checked easily from the expression ( |5.4| ). Notice that for < v < 1/2 
(1/2 < v < 1) the D2 brane passes through the point r = 0. This behaviour is similar to the 
p < 6 case. However, the present solution differs substantially from those studied in section 
4. In order to make this difference manifest, let us introduce the cylindrical coordinates 
(z, p) as follows: 



-r cos 9 



P 



r sin 9 . 



(5.6) 



The behaviour of the coordinate p in the region in which r — > oo is of special interest to 
interpret the asymptotic behaviour of the brane. It is not difficult to find how p depends on 
9 for the solution (15. 41): 



P 



* ( 9 
A tan - 

v 2 



i 

l-2i/ 



(5.7) 



Clearly, when u < 1/2, p ^ oo for 9 — » tt, whereas if v > 1/2 the coordinate p diverges for 
9 w 0. This implies that our solution does not behave as a tube |. This fact is more explicit 
if we rewrite our solution ( |5.4| ) as a function z = z(p). After a short calculation, we get: 



P 



2u 



2 A 1 - 2 " 

When p — > oo, the function z(p) behaves as: 



(p2(l-2^) _ A 2(l-2u) ^ 



(51 



P 



2(1-") 



if ^ < 5 , 



lim zip 



-p 



2// 



if v > i 



(5.9) 



2 ' 



5 A similar calculation for p < 6 gives p — > constant at the spikes. 
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(a) (b) 

Figure 3: Representation of the D2 embedding for v = 1/4 (a) and v — (b). The conven- 
tions are the same as in figure 2. 



and thus, as expected, z — > +oo (z — > — oo) for z/ < 1/2 (z/ > 1/2). Moreover, eq. Q5.9Q 
implies that the asymptotic shape of the solution is that of a paraboloid Q rather than a tube. 
We have thus an "upper paraboloid" for v < 1/2 and a "lower paraboloid" for v > 1/2. In 
figure 3 we have plotted the p = 6 solution for different values of v. Notice that for v ^ 0,1 
z(p) has an extremum for some value of p ^ 0, while for u — 0,1 the extremum is located 
at p = 0. This extremum is a minimum (maximum) for v < 1/2 [y > 1/2). It is interesting 
to point out that, although the solution for p = 6 differs from the ones found for p < 5, the 
energy U is still given by \Z S \, where Z s is the same as in eq. ( f4.30|) . 



It follows from our previous analysis that v = 1/2 is a critical point in the behaviour 
of the D6-D2 system. The BPS condition ( |5.2|) is not defined in this case and we have to 
come back to the equation of motion (|3.1|) . Fortunately, we will be able to find the general 
solution of the field equation in this case. The key point in this respect is the observation 
that, for p = 6 and v = 1/2, the near-horizon energy density does not depend on 9 explicitly. 
Indeed, a simple calculation shows that: 

U = T 2 Q 1 R J d9Vr 2 + r' 2 , („ = -). (5.10) 

This ^-independence implies the "conservation law": 

dU 

r' — — — U = constant , (5.11) 
or ' 



or, using the explicit form of U given in eq. (5.1U): 



4 Strictly speaking, only for v = 0, 1 eq. (|5.8|) represents a paraboloid 
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I 
I 



Figure 4: Solution of the equation of motion for p = 6 and v — 1/2 for the near-horizon 
metric (eq. ( |5.13| )). 



,j. 2 



= = C , (5.12) 

where C > is constant. It is not difficult to integrate the first-order equation ( |5.12| ). The 
result is: 

r(0) = ■ 7^7 S , (5-13) 

a being a new integration constant. The solution (|5.13|) represents a cone with opening 
half-angle a and with its vertex located at a distance C/sma from the origin (see figure 4). 
This fact is neatly shown if we rewrite ( 5.13Q in terms of the cylindrical coordinates (p, z) 



defined in eq. ( |5.6| ). In these coordinates, eq. ( |5.13| ) becomes 



pcosa — z since = C . (5-14) 

It is clear from figure 4 that a can take values in the range —tc/2 < a < it/2. When a = 
the cone degenerates in the cylinder p = C, while for a = ±7r/2 the solution (|5.14 ) is the 



plane z = =pC. Notice that eq. Q5.14 ) is invariant under the transformation a — > —a, 



z — > —z. Thus, for a > (a < 0) z goes to +oo (— oo) when the angle 6 takes the value 
7i — a (—ex). It is interesting to point out that, when p — > oo, z diverges linearly with p, 
which is, certainly, the limiting case of the v ^ 1/2 situation (see eq. ( |5.9| )). This softer 
behaviour for v = 1/2 will be also reflected in the solution for the full Dp-brane metric, as 
we shall check in the next section, where we will find the solutions of the BPS equation for 
the full Dp-brane brane metric ( j2.1|) . We will check that for p = 6 and v — 1/2 this solution 
is equivalent, in the near-horizon region, to our general solution Q5.14 ) with C — 0, i.e. only 



those cones passing through the origin correspond to BPS solutions of the asymptotically 
flat metric. 
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6 Asymptotically flat background 



In this section we are going to study the solutions of the BPS condition beyond the near- 
horizon approximation. Recall that the full metric of the Dp-brane background is obtained by 
taking a = 1 in the harmonic function ( |2.2|) . The full Dp-brane metric is asymptotically flat 
and, therefore, one expects that a brane placed far away from the location of the background 
branes is not bent by the gravitational field. We shall concentrate first in the study of this 
asymptotic behaviour and, afterwards, we will consider the complete solution. 

It is convenient for our purposes to rewrite the BPS condition ( |3. 7[ ) in the coordinates 
(z,p) of eq. (151)1) . m these coordinates, the D(8-p)-brane embedding is characterized by a 
function z(p). The BPS equation gives dz/dp in terms of (z,p). Actually, it is not difficult 
to relate z(p) and its derivative to r(9) and r'{6). By using only the relation between both 
coordinate systems one can verify that: 

dz sin 9 - cos 9 ^ 

— = — ■ o.l 

dp cos9 + sm9 ^ v ; 

By substituting the value of the ratio r'/r given by eq. (|3.7p for a = 1 on the right-hand side 
of this equation, one arrives at the following first-order differential equation for the function 

z(p): 



dp " B7-P + [p 2 + z 2 ] 1 ^ P 7 - p 



dz \p 2 + z 2 }^ E D p [aicttw(-p/z, j 

(6.2) 

It is now elementary to evaluate the right-hand side of eq. ( |6.2j ) in the asymptotic region in 
which p — > oo, z/p — > (and 9 — > tt/2). In this limit, eq. (|6.2|) takes the form: 

* (6 .3) 
dp p' p 

where we have only kept the first term in the expansion in powers of R/p. Notice that 
z'(p) — > as p — > oo, according to our expectations. However, also in this approach, the 
behaviour of the p = 6 case differs from that corresponding to p < 6. Indeed, for p — 6, 
z'(p) ~ and z(p) diverges logarithmically in the asymptotic region while, for p < 6, z(p) 
approaches a constant value as p — > oo. Let us specify further these two kinds of behaviours. 
The coefficient multiplying the power of p on the right-hand side of eq. ( |6.3p is D p (ir/2), 
which can be computed from eq. (|2.22|) : 



D p (n/2) = -2^ %^Jr R 7 - p ( I - v ) . (6.4) 



r(V) 2 

Using this result we can integrate eq. ( |3.3|) . For p < 6, we get: 
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where is a constant representing the asymptotic value of z. From eq. ( |6.5| ) one concludes 
that the sign of z^ — z depends on the sign of ~ — v. i.e. if v < 1/2 (y > 1/2) the 
brane reaches its asymptotic value of z from below (above). If p = 6 we get the expected 
logarithmic dependence: 

z(p) ~ (1 - 2u)R\ogp + ■•• . (6.6) 

Notice that when p — * oo in eq. ( |6.6|) , -2 — > +oo if i/ < 1/2 while, on the contrary, z — > — oo 
for f > 1/2, in agreement with our near-horizon analysis of section 5. It is interesting to 
point out that, although in this case z diverges as p — > oo, z/p still vanishes in this limit, as 
assumed in the derivation of eq. ( |6.3| ) . 

For v = 1/2, the leading term in the asymptotic expansion vanishes and we have to 
compute the next-to-leading term. This does not change significantly the analysis of the 
p < 6 case. However, for p = 6, things change drastically when v — 1/2. In fact, one can 
prove from eq. ( |6.2| ) that, in this case, z approaches a constant value z^ as p — > oo (and, 
actually, z — z^ decreases as for p » R). 

Let us now see how one can integrate exactly the BPS condition for the full metric. 
The resulting solution must reproduce the asymptotic behaviour we have just described. 



For p = 3, eq. (|6.2|) was integrated numerically in ref. |Tl|]. As we shall show below, our 
analytical solution agrees with these numerical results. It is more convenient to come back 
to our original (r, 9) coordinates. The basic strategy we will adopt to integrate the BPS 
equation is to write it in terms of the same function A p ( 9 ) that we have already used to find 
the near-horizon solutions. For illustrative purposes, we will do it for a general value of the 
parameter a of the harmonic function. After using the definition ( [Op , eq. fl3.7| ) can be put 
in the form: 

r ' _ a r 7 " p ( sin 9 ) 8 ~ p + R 7 ~ p ( sin 9 f' p - R 7 ~ p A p ( 9 ) cos 9 

7 ~ a r 7 ~ p ( sin 9 ) 7 ~ p cos 9 + R 7 ~ p A p ( 9)sm9 ' ^ ' ' 



Notice that eq. (|6.7|) for a = is identical to eq. (f4.2|) . Let us now isolate the terms of (|6.7| ) 
depending on a in one of the sides of the equation. Doing this, it is elementary to prove that 
eq. (|6.7p can be rewritten as: 



ar 7 - p (sm9) 7 - p ^- (rcostf) = R 7 ~ p \r ( sintf f~ p - A P (9)^- (rsin^)l . (6.8) 
d9 L 

On the right-hand side of this expression we recognize a term containing (sin6') 6_p , whose 
representation as a derivative of the function A p ( 9 ) was crucial to integrate the near- horizon 
BPS condition for p ^ 6. Actually, assuming that p ^ 6 and using eq. (|4.4j) , our differential 
equation becomes: 

ar 7 ~ p (sm9) 7 ~ p ^- (rcostf) = 
d9 

= I*lJL [ rsin e L A (6-p)A p (9) (rsm9)} , (p^6) . 

6 — p 1 a9 d9 J 

(6.9) 
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v = 



v = l 

2 



Figure 5: Solutions of the BPS differential equation for the asymptotically flat metric 
(eq. ( |6.12j )) for p = 3 and several values of v and z^. 



Taking a = we recover the solutions ( [4.5| ) and ( f4.18| ). From now on we shall put a = 1 in 
all our expressions. In this case we can pass the r 7 ~ p ( sin# ) 7_p factor in the left-hand side 
of eq. (|6.9|) to the right-hand side and, remarkably, this equation can be written as: 



d 
d9 



r cos 9 



R 7 ~p d 



6 — p d6 



A p (9] 



r sin 9 



\p^6 



The integration of eq. (|6.10|) is now immediate: 



r cos# 



R 7 ~p 



A p (e; 



6 — p ( r sin 9 ) 6 p 



- z n 



p^6 



(6.10) 



(6.11) 



where is a constant of integration which can be identified with the asymptotic value of z 
introduced in eq. (|6.5|). Actually, in terms of the coordinates (z,p), our solution ( |6.11| ) can 
be written as: 



_r 7 ~p Ap ( arctan (—p/z) 



'p^6 



(6.12) 



6 - p p b -P 

Eqs. ( p. 11| ) and (|6.12 ) give the analytical solution of the BPS differential equation for p ^ 6 

oo is reduced to eq. 



It is easy to verify that eq. (|6.12|) for p — » oo is reduced to eq. (|6.5|) . Moreover, by derivating 
both sides of eq. ( |6.12| ) with respect to p, one can extract the value of dz/dp. The result in 
terms of (z, p) is precisely the one written in Q6.2p , i.e. the function z(p) defined by (|6.12|) is 
a solution of the BPS condition in cylindrical coordinates. 

Let us now study the analytical solution we have found and, specially, how it interpolates 
between the near-horizon and asymptotic regions. Notice, first of all, that eqs. ( |6.11| ) and 
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( 3.12j ) give r(6) and z(p) in implicit form. In particular, eq. ( 3.11| ) determines r(6) as the 
root of a polynomial of degree 7 — p which, for p > 3, can be solved algebraically. However, 
the implicit equations ( |6.11| ) and ( |6.12| ) are much easier to handle and, actually, can be 
used to plot the embeddings. In figure 5 we have represented these embeddings for different 
values of v and z^. Moreover, by studying eqs. ( |6.11| ) and ( |6.12| ) we can characterize the 
general features of our BPS solutions. The detailed analysis of these solutions is performed 
in appendix A. It follows from this analysis that our solution ( |6.12[ ) coincides with the near- 
horizon embeddings in a region close to the origin. Actually, as can be seen from the plots of 
figure 5, the integral of the BPS differential equation for the full metric contains cylindrical 
regions which resemble closely the tubes found in the near-horizon analysis. It is not difficult 
to relate the parameter z^ of eq. ( |6.12|) to the constant C appearing in the lower and upper 
tube solutions (eqs. ( |4.5|) and ( |4.18|) ). For large | ^ | it can be seen that both types of 
solutions are approximately equal if we are sufficiently near the origin and if the following 
relation: 



R 7 ~p 



(6.13) 



( 6 - p ) | I 
is satisfied. 

For < v < 1 the solution reaches the origin r = by means of a tube, which is an upper 
or lower tube depending on the sign of (see figure 5). For | z^ \ — > oo we will verify that 
this tubes can be regarded as bundles of fundamental strings, exactly in the same way as in 
the near-horizon case. 

The behaviour of the v = solution is different from that we have just described for 
< v < 10. When z^ — + — oo, the D(8-p)-brane is nearly flat. As we increase z^, an upper 
tube develops and for z^ — > +oo a bundle of N fundamental strings, connecting the Dp 
and D(8-p) branes, is created. As pointed out in ref. |Ll]], this solution provides a concrete 
realization of the so-called Hanany-Witten effect [[HJ. 

Let us now study the energy of our solutions. It is clear from the plots of figure 5 
that r is not a single-valued function of 9. It is thus more convenient to use the (z, p) 
coordinates in order to have a global description of the energy functional. From eq. ( |2.18| ) 
it is straightforward to obtain the expression of U as an integral over p. The result can be 
put in the form: 



U 



7-p 



dp 



A, 



D - 

p dp 



+ D 



. dz 

A 'T P 



(6.14) 



and A p is the function defined in 



In eq. ( |3.14j ) D p depends on p and z(p) as in eq. ( |6^ 
eq. d3.2| ). It is important to stress that eq. ( |6.14| ) gives the energy for any embedding z(p). 
Notice that on the right-hand side of eq. ( |6.14| ) we have the sum of two squares. The second 
of these two terms vanishes when the BPS condition dz/dp 
It is clear from these considerations that, if we define: 



-D p /A p (see eq. (gj)) holds 



X EE A r 



D, 



dz 
dp 



(6.15) 



3 The solution for v = 1 can be reduced to that for v = (see eq. ( A.l )) 
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we obtain the following lower bound for the energy of any embedding: 



U > T 8 _pfi 7 „ p J dp 



X 



(6.16) 



The bound ( |6.16D is saturated precisely when the embedding satisfies the first-order BPS 
equation Q6.2|) . A remarkable aspect of the function X appearing in the right-hand side of 
eq. (|6.16p is that it can be put as a total derivative. Indeed, one can verify that: 



X = — 



d 
dp 



-zD r 



R 



7-p 



V 



p 2 + z 1 



2 



,8-p 



(6.17) 



Eq. (|6.17|) is the analog in these coordinates of eq. (|3.8|) . As it happened with eq. (|3] 



in order to prove eq. (|6.17| ) one only has to use eq. (|2.14|) and, therefore, ( |6.17p is valid for 
any function z(p). An important consequence of eq. (|6.17 ) is that the bound (|6.16|) only 
depends on the boundary conditions of the embedding at infinity. Thus, one can say that 
the BPS embeddings are those that minimize the energy for a given value of z(p) at p — ► oo. 
From eq. ( |6.14[ ) it is very easy to obtain the energy Ubps of one of such BPS solutions. The 
result one arrives at is: 



U 



BPS 



dp 



dz / dz \ -i 

- + (t p ) 



dp 



D r 



(6.18) 



Apart from being simple, eq. ( |6.18| ) is specially suited for our purposes. Notice, first of all, 
that the integrand on the right-hand side of (|6.18|) is always non-negative due to eq. (|6.2|) . 
Secondly, we can use ( |6.18| ) to evaluate the energy of a tubular portion of the brane. Indeed, 
for large |zoo|, ^ is large in one of these tubes and the argument of D p is almost constant 

and equal to tt (0) for an upper (lower) tube. Neglecting the term containing ( ^ ) on the 
right-hand side of eq. ( |6.18p , and defining the length of the tube as: 



■'tube 



tube 



dp 



dz 



dp 

we get the following values for the energy of the tubes: 



(6.19) 



- T 8 _ p tt 7 _ p D p (tt) L tube , (upper tube) , 

r,„,„ = ! (6.20) 

8-p &7-p D p ( ) L tube , (lower tube) . 

Taking into account eqs. (|2.23 ) and (12.25 ), and using eq. Q2.24T ), one can easily show 
from the result in eq. ( |6.20| ) that Utube coincides with \Z S \, where the values of Z s are given 
in eq. (|4.30| ). This confirms our interpretation of the tubes as bundles of (1 — u)N and uN 
fundamental strings. 

It is not difficult to calculate the energy of the whole brane. Actually, by means of 
eq. ( |6.17| ) we can perform the integral appearing in the right-hand side of eq. fl6.18|) . As we 
are calculating the energy of an infinite brane, this integral is divergent. In order to regulate 
this divergence, let us introduce a cutoff p c , in such a way that the integral in eq. ( |6.18| ) is 
performed between p = and p = p c . It is not difficult to check from the properties of our 
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solutions that the contribution of the lower limit p = is zero and, thus, only the value at 
the cutoff contributes. After this process, Ubps can be put as: 



U. 



BPS 



—T 8 _ p VLj_ p 



zD r 



V 



+ 



R 7 ~p 



p 2 + z 2 



P 



8-p 



(6.21) 



P = Pc 



We have now to take p c — > oo and, therefore, z(p c ) — > z^. Notice that, in this limit, the 
argument of the function D p is ir/2. Using eq. (|6.4j ), Ubps can be put as a sum of two terms, 
namely: 



u 



BPS 



- - u)NT f z. 



Pc 



+ R 7 ~ p p, 



V 



(6.22) 



Notice that the first term on the right-hand side of eq. ( |6.22| ) is finite and depends linearly 
on Zoo, while the second term diverges when p c — > oo and is independent of z^. According 



to ref. |TB[, one can give the following interpretation to this divergence. Let us consider a 
D(8-p)-brane embedded in the metric fl2.1|) along the plane z = or, equivalently, such that 
its worldvolume is determined by the equation 6 = n/2. Notice that in this configuration 
the brane is not bent at all and, for this reason, it will be referred to as the "ground state" 
of the brane. Let g gs be the induced metric on the worldvolume of the D(8-p)-brane for this 
ground state configuration. Putting the worldvolume gauge fields to zero and substituting 
g by g gs in eq. (|2.6|) , we get the action of the ground state: 



S, 



8-p 



-det (g gs ) 



The energy E gs of the ground state is obtained from S gs as: 



E n 



S, 



gg 
T 



(6.23) 



(6.24) 



where T = J dt. Using the metric given in eq. ( |2.1| ) and the value of the dilaton field 
displayed in eq. (2.5), the calculation of E gs is a simple exercise. By comparing the result of 



this computation with the right-hand side of eq. ( p. 22] ), one discovers that E gs is equal to 
the divergent contribution to Ubps- Therefore, one can subtract E gs from Ubps an d define 
a renormalized energy U ren as: 



U r . 



u, 



- E n 



'BPS ^gs 

It follows from eq. (|6.22|) that U ren is given by: 



NTf z c 



(6.25) 



(6.26) 



By derivating U ren with respect to z^, we learn that there is a net constant force acting on 
the D(8-p)-brane, which is equal to ( | — u^j NTf, i.e. equivalent to the tension of ( | — v J Af 
fundamental strings. In ref. [|llj this force was interpreted as a consequence of the fact that, 
due to the p + 2-form flux captured by the D(8-p)-brane, the latter is endowed with an 
effective charge equal to ( ~ — v ) N units. 



23 



Figure 6: BPS embedding of the D2-brane in the full D6-brane geometry for v = 1/2. 



Let us now integrate the BPS differential equation for p = 6. Using the value of A 6 (0) in 
eq. Q6.81 ), we get (for a = 1): 



r sin 9 — ( r cos 9) = R 
d9 



r + ( 2v — 1 ) — ( r sin 9 
d9 



p = 6 



This equation can be recast in the form: 



(6.27) 



d_ 
dl 



r cos 9 



R 



dl 



log 



tan 



r sin 9 ) 



2u-l 



P 



6 



(6.28) 



and, therefore, its integration is immediate. It is more interesting to write the result in 
cylindrical coordinates. After a short calculation one gets: 



(z + ^p 2 + z 2 ) 



P 



P 



6 



(6.29) 

where A; is a constant of integration. It is easy to check that the function z(p) parametrized 
by eq. ( |6.29D has the asymptotic behaviour displayed in eq. (|6.6|) . It follows that, as 
expected, z(p) does not reach a constant value for p = 6 and v ^ 1/2. The profiles of the 
v 7^ 1/2 curves are similar to the near- horizon ones, plotted in figure 3. The only difference 
is that, instead of the behaviour written in eq. Q5.9p, | z(p) \ grows logarithmically as p — > oo. 
Moreover, similarly to what happens to the p < 5 solution, eq. ( |6.29|) is invariant under the 
t r ansf or mat ion : 



P ^ P- 



(6.30) 



It is also easy to check that the function z(p) defined by eq. ( |6.29| ) coincides, in the 
near-horizon region, with the one written down in eq. ( |5.8|) if the constants A and k are 
identified as follows: 



A 



2u-l 



1 



(6.31) 



For v — 1/2 the behaviour of the solution changes drastically. Indeed, in this case, one 
can solve eq. ( |5.29| ) for p and get p = p(z). The result is: 
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P = — k _ z , (u = -,p = 6). (6.32) 



sinh^ ' 2 



The function ( |6.32|) has been represented graphically in figure 6. From this plot one notices 
that linip^oo z(p) = k. Actually, one can prove using eq. (|6.32j ) that for k > (k < 0) the 
variable z takes values in the range 0<z<k(k<z<0). Thus, in this (p = 6, v = 1/2) 
case, the constant k plays the same role as z^ did for p < 5. In particular, for k = the 



solution (|6.32|) is equivalent to the equation z = 0, which is precisely the same solution we 



found for p < 5, v = 1/2 and z^ = 0, i.e. one obtains the "ground state" solution in both 
cases. 

Let us finally point out that, by linearizing eq. (|6.32|) in p and z, one gets the equation 



of a cone passing through the origin and with half-angle equal to: 



1 

tana = -t— -y . (6.33) 

sinh (J) 



Therefore, with the identification ( |6.33D , we get a perfect agreement with the near- horizon 
solutions fl5T|) with C = 0. 



7 Summary and Conclusions 

In this paper we have studied the embedding of a D(8-p)-brane in the background geometry 
of a stack of coincident Dp-branes. This embedding is governed by the worldvolume action of 
the D(8-p)-brane (eq. fl2.1|) ), which determines the equation of motion (|3.1| ). By using a BPS 
argument we have found a bound for the energy of the system such that those embeddings 
which saturate it are also a solution of the equation of motion. This equation of motion 
(eq. fl3~7f) ) is a first-order differential equation which, amazingly, can be solved analytically 
both in the near-horizon and asymptotically flat geometries. 

The solutions of the BPS equations give the deformation of the D(8-p)-brane under the 
influence of the gravitational and RR fields created by the background branes. Generically, 
these solutions contain tubes connecting the D(8-p)-brane to the background branes which, 
after analyzing its energy, can be interpreted as bundles of fundamental strings. From the 
point of view of the gauge theory defined on the worldvolume of the background branes, 
these tubes represent baryonic multiquark states. 

There are several topics which were not covered by our analysis which, in our opinion, 
it would be interesting to consider in a future work. Let us mention some of them. First 
of all, it would be interesting to understand the relation between the BPS condition and 



super symmetry. One expects, following the line of ref. |C|, that the BPS differential equation 
is precisely the requirement one must impose to the brane embedding in order to preserve 
some fraction of the space-time supersymmetry. An analysis of the supersymmetry algebra, 



by applying the methods of ref. pi |, could shed light on this aspect of the problem. 

Another topic which would be worth to explore is the use that our exact results could have 
in the study of baryons in gauge theories both in the supersymmetric and non-supersymmetric 
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models. In the latter case we have to extend our results to the situation in which the back- 
ground brane configuration is not extremal. 

Finally, it would be interesting to find out if there exist BPS conditions, similar to the 
ones found here, for brane embeddings in more general brane geometries, such as the ones 
corresponding to the intersection of several branes of different types. 

The final goal of these studies is twofold. On the one hand, we would like to know what 
string theory can teach us about the non-perturbative structure of gauge theories while, at 
the same time, we would like to uncover aspects of a more complete formulation of string 
theory. 
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APPENDIX A 

In this appendix we are going to carry out an study of the solution found in section 6 of 
the BPS differential equation for the asymptotically flat metric for p < 5 (eqs. ( |6.11| ) and 

(ill))- 



The first thing we want to point out in this respect is that, due to the property of A p 
displayed in eq. (|4.19|) , our solution ( 6.12| ) is invariant under the transformation: 



p^p, z-^-z, z^-y-z^, v ► 1 v . (A.l) 

Due to this invariance we can restrict ourselves to the case < v < 1/2. The solutions 
outside this range of v can be obtained by performing a reflection with respect to the z = 
axis. Therefore, unless otherwise stated, we will assume in what follows that v < 1/2. 
Moreover, it is more convenient for our purposes to rewrite eq. (|6.12 ) in the equivalent form: 



R 7 ~ p / 

z ) p 6 ~ p = A p ( arctan ( -p/z) ) . (A.2) 



6 — p 

We will start our analysis of the solution by studying its cut with the p = axis. By looking 
at eq. (|A.2[) it follows immediately that A p must vanish for p = 0. Thus, the angle 9 at which 



the solution reaches the p = axis must be 9 = 6q, 9q being the same angle appearing in 
the near-horizon solution (see eq. (|4.8|) ). Following eq. (|5.6j ), p = r(6) sin 9 and, therefore, 
if p = for 9 = 9q one must have r(9o) sin#o = 0. There are two possibilities to fulfill this 
equation. If, first of all, v ^ 0, the angle 9 is non-vanishing, which means that, necessarily, 
r(6 l ) = 0, i.e. the solution reaches the origin r = at an angle 9 . Notice that, as it should 
occur near r = 0, this is precisely what happens in the near-horizon solution. The second 
possibility is v = 0. In this case 9o = and the vanishing of p does not require r = 0. Again, 
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this is in agreement with the near-horizon solution. Actually, the distance r(0) at which the 
v = solution cuts the p = axis can be determined, as in section 4, by expanding A p (8) 
around 9 = 0. Using eq. ( }4.15| ) one gets that r(0) must satisfy: 

[r(0)] 6 " p [r(0) + Zoo ] = |— , {" = 0). (A.3) 

o — p 

As the right-hand side of this equation is positive, it follows that r(0) > — z^. When this 
condition is satisfied, the left-hand side of eq. (|A.3| ) is a monotonically increasing function 
of r(0) and, therefore, there exists a unique solution for r(0). 

Let us now consider the cut with the z = axis. We shall denote the corresponding value 
of the p coordinate by p , i.e. po = p( z — 0). If p 7^ 0, the value of the angle 9 for z = 
is 9 = Ti / 2 and, after evaluating the function A p for this value of 9, eq. ( |A.2j ) gives: 

Notice that for v ^ 1/2 the right-hand side of eq. ( |A.4|) only makes sense for z^ > 0. 
Therefore, only for > the solution for v ^ 1/2 cuts the 2 = axis at a finite non- 
vanishing value of the coordinate p. For u = 1/2 and ^ ^ eq. ( [A.4] ) has no solution for 
Po 7^ 0. All these features appear in the plots of figure 5. 

From eq. QA.2 ) one can extract the range of allowed values of the coordinate z. Indeed, it 



follows from eq. (|A.2 ) that the signs of z^ — z and A p are the same. On the other hand, we 



know that A p (6) is positive for 9 > 9 and negative for 9 < 9 . Thus, when 6 > 6 (6 < 6 ) 
one must have z < z^ (z > z^), while for 9 = 9 either p = or else z = z^. As 6 < n/2 
for v < 1/2, one can easily see that these results imply that z < for z^ < 0, whereas z can 
be positive or negative if z^ > 0. 

It is easy to prove that z, as a function of p, must have a unique extremum. For v < 1/2 
this extremum is actually a minimum, as we are going to verify soon. From eq. ( |6.2j ) it is 
clear that dz/dp vanishes if and only if D p is zero. Recall that the Gauss' law (eq. (|2.14| )) 
implies that D p is a monotonically decreasing function of 9 (see figure 1). Moreover, it follows 
from eqs. ( |2.23| ) and ( |2.25|) that D p (0) > and D p (tt) < 0. Thus, it must necessarily exist 
a unique value 9 m of 9 such that: 

D p (9 m ) = 0. (A.5) 

Clearly, at the point of the curve z(p) at which 9 = 9 m the derivative dz/dp is zero. It is 
interesting to compare 9 m with the angle 9 for which A p is zero. By substituting 9 = 9$ va. 
the equation which relates D p ( 9 ) and A p ( 9 ) (eq. (fOT)), one gets: 

D p ( 9 ) = R 7 ~ p ( sin 9 f~ p cos 9 . (A.6) 

For v < 1/2 one has 9 < n/2 and, thus, eq. ( |A.6| ) gives D p (9 ) > 0. The monotonic 
character of D p (9) implies that < 9 < 9 m < n/2. Notice that 9q = 9 m if v = 
($o — dm = 0) or v = 1/2 (#o = 9 m = vr/2). The coordinate z of the extremum is 
z m = — r (#m) co$9 m . As 9 m < tt/2, one must have z m < 0. Moreover, since 9 m > 9 Q , then 
A p (9 m ) > and eq. ( |A.2| ) gives z m < z^. By using the value of A p at 9 = 9 m , which is: 
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A p (6 m ) = (sm6 m ) 6 - p cos# m , (A.7) 
(see eq. ( f4.1|) ) one can obtain an expression which determines z m , namely: 



1 6 — p 



R 7 ~ P ( „ \7~P 



,cos# m . (A.8) 
6 — p v ' 

It is not difficult to verify now that, when v < 1/2, the extremum at 9 = 8 m is a minimum. 
In order to prove it, one must evaluate d 2 z/dp 2 . This can be done by derivating eq. (|6.2| ). 



After putting p = p m = r(6 m ) sin 6 m , and using the fact that dz/dp vanishes for p = p m , 
one arrives at: 



d 2 z 



dp 2 



7 - v ) R 7 ~p 

' z m . (A.9) 



Pm Z m ) ^ P + ( Pm Z \ 



2 7 ~ 



, 2 
m 1 ~m I 



When v < 1/2, one has z m < and, as a consequence, the right-hand side of eq. ( |A.9| ) is 
strictly positive. Therefore z(p) has a minimum for p = p m as claimed. Notice that, when 
v = 1/2, i.e. d m = tt/2, the right-hand side of eq. ( |A.8| ) vanishes and, thus, z m = (the 
value z = Zoo is reached asymptotically). 



Eq. ( |A.8| ) can be solved immediately if z^ = 0. Indeed, in this case, one has: 



\7-p 



R 7 ~ P t „ \7 



6 — J9 



-p 



cos^ m , ( Zoo = 0). (A.10) 



For general values of z^ we cannot give the general solution of eq. (|A.8|) . However, one can 
extract valuable information from this equation. An important point concerning eq. ( |A.8| ) 
is the fact that its right-hand side does not depend on ^ (it only depends on p and v) and, 
therefore, it remains constant when we change the asymptotic value of the z coordinate. 
This allows us to study the behaviour of the minimum in different limiting situations. Let us 
suppose, first of all, that z^ — > +oo for u < 1/2. As z m < 0, the difference z^ — z m — » +oo. 
In order to keep the right-hand side of eq. (|A.8| ) finite, one must have z m — > 0~. Actually 
z m should behave as: 

1 1 

( Zoo -»• +oo, v < -) . (A.ll) 



( ^oo ) & - p 

On the other hand if z^ — > — oo, again for v < 1/2, the only possibility to maintain 
I z m | 6_p ( Zoo — z m ) constant is that — z m — > + (recall that z m < Zoo and, thus, z m — > 
is impossible if Zoo — > — oo). The actual behaviour of z^ — z m in this limit is: 

Zoo ~ Z m ~ l — — , ( Zoo -> -oo , 1/ < ^ ) . (A.12) 

\ Zoo ) ^ 

The v = embedding presents some characteristics which make it different from the 
v ^ ones. We have already shown that the v = solution cuts the p = axis at a 
coordinate z = — r(0), with r(0) determined by eq. ( |A.3|) . Moreover, for v = the function 



A p on the right-hand side of eq. (|A.2|) is non-negative for all p ^ 0, which implies that 



z < -Zoo for these solutions. Recall that, in this case, z(p) has a minimum at p = 0. For 
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Zoo — > —oo, due to eq. ( |A.12| ), z(p) is approximately a constant function, i.e. z(p) ~ z c 
If, on the contrary, z^ — > +00, an "upper tube" connecting the point z = — r(0) — > (see 
eq. ( | A . 1 1| ) ) and the asymptotic region is developed (see figure 5). 



v = 1/2 is also a special case. Indeed, for v = 1/2 the angle 60 is equal to ir/2 and, 
therefore, the function A p in eq. ( |A.2| ) takes positive (negative) values for z > [z < 0). By 
studying the signs of both sides of eq. (|A.2j ) one easily concludes that for > one must 



have < z < Zoo, whereas for z^ < the coordinate z takes values on the range z^ < z < 
(see figure 5). When z^ = 0, the analysis of eq. (|A.2 ) leads to the conclusion that the 
solution is the hyperplane z = 0. 

It is also possible to determine the range of possible values that the coordinate 9 can take. 
Actually 9 has generically an extremal value, which we shall denote by 9*. The existence of 
9* can be proved by computing the derivative d9/dr for our solution. A short calculation 
proves that this derivative vanishes if the coordinate z takes the value: 

z* = — ^oo • (A. 13) 

7 — p 

Alternatively, one could determine z* as the point at which the denominator of the BPS 
condition ( |3.7| ) vanishes. By computing the second derivative d 2 9/dr 2 at 9 = 9*, one can 
check that 9* is a minimum (maximum) if z^ < (zoo > 0)Q. It is also interesting to point 
out that for z = z* the function Z of eq. Q4.21Q changes its sign. This is consistent with the 
extremal nature of 9*. 
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